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Abstract 

We propose a formulation of the term structure of interest rates in which the forward 
curve is seen as the deformation of a string. We derive the general condition that the 
partial differential equations governing the motion of such string must obey in order to 
account for the condition of absence of arbitrage opportunities. This condition takes a 
form similar to a fluctuation-dissipation theorem, albeit on the same quantity (the for- 
ward rate), linking the bias to the covariance of variation fluctuations. We provide the 
general structure of the models that obey this constraint in the framework of stochastic 
partial (possibly non- linear) differential equations. We derive the general solution for 
the pricing and hedging of interest rate derivatives within this framework, albeit for 
the linear case (we also provide in the appendix a simple and intuitive derivation of 
the standard European option problem). We also show how the "string" formulation 
simplifies into a standard N- factor model under a Galerkin approximation. 



Short title : Forward Interest Rate dynamics 



PACS: 02.50.-r Probability theory, stochastic processes, and statistics 
05.40. +j Fluctuation phenomena, random processes, and Brownian motion 
89.90. +n Other areas of general interest to physicists 



1 Introduction 



Imagine that Julie wants to invest $1 for two years Q. She can devise two possible 
strategies. The first one is to put the money in a one- year bond at an interest rate r%. 
At the end of the year, she must take her money and find another one-year bond, with 
interest rate r\ which is the interest rate in one year on a loan maturing in two years. 
The final payoff of this strategy is simply (1 + + r\). The problem is that Julie 
cannot know for sure what will be the one-period interest rate r\ of next year. Thus, 
she can only estimate a return by guessing the expectation of r\. 

Instead of making two separate investments of one year each, Julie could invest her 
money today in a bond that pays off in two years with interest rate r^. The final payoff 
is then (1 + r^) 2 . This second strategy is riskless as she knows for sure her return. 
Now, this strategy can be reinterpreted along the line of the first strategy as follows. 
It consists in investing for one year at the rate t\ and for the second year at a forward 
rate f2- The forward rate is like the r\ rate, with the essential difference that it is 
guaranteed : by buying the two-year bond, Julie can "lock in" an interest rate fi for 
the second year. 

This simple example illustrates that the set of all possible bonds traded on the 
market is equivalent to the so-called forward rate curve. The forward rate f(t,x) is 
thus the interest rate that can be contracted at time t for instantaneously riskless 
borrowing Q or lending at time t + x. It is thus a function or curve of the time- 
to-maturity x 0, where x plays the role of a "length" variable, that deforms with 
time t. Its knowledge is completely equivalent to the set of bond prices P(t, x) at 
time t that expire at time t + x (see eq.(|j) below). The shape of the forward rate 
curve f(t, x) incessantly fluctuates as a function of time t. These fluctuations are due 
to a combination of factors, including future expectation of the short-term interest 
rates, liquidity preferences, market segmentation and trading. It is obvious that the 
forward rate f(t,x + Sx) for Sx small can not be very different from f(t,x). It is thus 
tempting to see f(t,x) as a "string" characterized by a kind of tension which prevents 
too large local deformations that would not be financially acceptable. This superficial 
analogy is in the follow up of the repetitious intersections between finance and physics, 
starting with Bachelier § who solved the diffusion equation of Brownian motion as a 
model of stock market price fluctuations five years before Einstein, continuing with the 
discovery of the relevance of Levy laws for cotton price fluctuations by Mandelbrot @ 
that can be compared with the present interest of such power laws for the description 
of physical and natural phenomena [||. We could go on and cite many other examples. 
We investigate how to formalize mathematically this analogy between the forward rate 
curve and a string. In this goal, we formulate the term structure of interest rates as the 
solution of a stochastic partial differential equation (SPDE) ||, following the physical 
analogy of a continuous curve (string) whose shape moves stochastically through time. 

The equation of motion of macroscopic physical strings is derived from conserva- 
tion laws. The fundamental equations of motion of microscopic strings formulated to 
describe the fundamental particles || derive from global symmetry principles and du- 

1 "Instantaneous riskless" describes the fact that the forward rate is the rate that applies for a small time 
increment 5t as seen from equation (||) below and is fixed during this time, thus being locally riskless. 

2 The maturity of a financial product is simply its lifetime. In other words, it is the time interval between 
the present and the time of extinction of the rights attached to the financial product. 
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alities between long-range and short-range descriptions. Are there similar principles 
that can guide the determination of the equations of motion of the more down-to-earth 
financial forward rate "strings"? 

The situation is a priori much more difficult than in physics as illustrated by the 
following picturial analogy quoted from the journalist N. Dunbar at Financial Products 
magazine. Suppose that in the middle ages, before Copernicus and Galileo, the Earth 
really was stationary at the centre of the universe, and only began moving later on. 
Imagine that during the nineteenth century, when everyone believed classical physics 
to be true, that it really was true, and quantum phenomena were non-existent. These 
are not philosophical musings, but an attempt to portray how physics might look if 
it actually behaved like the financial markets. Indeed, the financial world is such 
that any insight is almost immediately used to trade for a profit. As the insight 
spreads among traders, the "universe" changes accordingly. As G. Soros has pointed 
out, market players are "actors observing their own deeds". As E. Derman, head of 
quantitative strategies at Goldman Sachs, puts it, in physics you are playing against 
God, who does not change his mind very often. In finance, you are playing against 
Gods creatures, whose feelings are ephemeral, at best unstable, and the news on which 
they are based keep streaming in. Value clearly derives from human beings, while mass, 
charge and electromagnetism apparently do not. This has led to suggestions that a 
fruitful framework to study finance and economy is to use evolutionary models inspired 
from biology and genetics. 

This does not however guide us much for the determination of "fundamental" equa- 
tions, if any. Here, we propose to use the condition of absence of arbitrage opportunity 
and show that this leads to strong constraints on the structure of the governing equa- 
tions The basic idea is that, if there are arbitrage opportunities (free lunches), they 
cannot live long or must be quite subtle, otherwise traders would act on them and 
arbitrage them away. The no-arbitrage condition is an idealization of a self-consistent 
dynamical state of the market resulting from the incessant actions of the traders (ar- 
bitragers). It is not the out-of- fashion equilibrium approximation sometimes described 
but rather embodies a very subtle cooperative organization of the market. 

We consider this condition as the fundamental backbone for the theory. The idea 
to impose this requirement is not new and is in fact the prerequisite of most models 
developed in the academic finance community. However, applying it in the present 
context is new. Modigliani and Miller 0, ^] have indeed emphasized the critical role 
played by arbitrage in determining the value of securities. It is sometimes suggested 
that transaction costs and other market imperfections make irrelevant the no-arbitrage 
condition Q. Let us address briefly this question before presenting our results. 

Transaction costs in option replication and other hedging activities ^ have been 
extensively investigated since they (or other market "imperfections") clearly disturb 
the risk-neutral argument and set option theory back a few decades. Transaction costs 
induce, for obvious reasons, dynamic incompleteness, thus preventing valuation as we 
know it since Black and Scholes ju]]. However, the most efficient dynamic hedgers 

3 Arbitrage, also known as the Law of One Price, states that two assets with identical attributes should 
sell for the same price and so should the same asset trading in two different markets. If the prices differ, a 
profitable opportunity arises to sell the asset where it is overpriced and to buy it where it is underpriced. 

4 Finance is all about risks but some risks can be hedged, i.e. offset, by trading in different financial 
instruments. 
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(market makers) incur essentially no transaction costs when owning options [] (see 
the appendix for a definition of options). These specialized market makers compete 
with each other to provide liquidity in option instruments, and maintain inventories 
in them. They rationally limit their dynamic replication to their residual exposure, 
not their global exposure. In addition, the fact that they do not hold options until 
maturity greatly reduces their costs of dynamic hedging. They have an incentive in the 
acceleration of financial intermediation. Furthermore, as options are rarely replicated 
until maturity, the expected transaction costs of the short options depend mostly on the 
dynamics of the order flow in the option markets -not on the direct costs of transacting. 
The conclusion is that transaction costs are a fraction of what has been assumed to be 
in the literature [11]. For the efficient operators (and those operators only), markets 
are more dynamically complete than anticipated. This is not true for a second category 
of traders, those who merely purchase or sell financial instruments that are subjected 
to dynamic hedging. They, accordingly, neither are equipped for dynamic hedging, nor 
have the need for it, thanks to the existence of specialized and more efficient market 
makers. The examination of their transaction costs in the event of their decision to 
dynamically replicate their options is of no true theoretical contribution. 

A second important point is that the existence of transaction costs should not be 
invoked as an excuse for disregarding the no-arbitrage condition but rather should be 
constructively invoked to study its impacts on the models. 

This work expands our previous work [12] which introduced stochastic strings mul- 
tiplied by volatility functions to shock forward rates. Our present approach directly 
deals with the SPDE equation for the forward rates in contrast to the shocks that 
drive the forward rate curve that were treated in Ref.[^]. Our model provides a fur- 
ther extension of the term structure model of Heath, Jarrow and Morton [14] and is as 
parsimonious and tractable as the traditional HJM model, but is capable of generating 
a much richer class of dynamics and shapes of the forward rate curve. Its main mo- 
tivation is to address the interplay between external factors represented by stochastic 
components (noise) and possible non- linear dynamics. 



2 Definitions 

We postulate the existence of a stochastic discount factor (SDF) that prices all assets 
in this economy and denote it by M. This process is also termed the pricing kernel, 
the pricing operator, or the state price density. We use these terms interchangeably. 
Ref.|15] is an excellent reference for the theory behind the SDF. It is well known that 



assuming that no dynamic arbitrage trading strategies can be implemented by trading 
in the financial securities issued in the economy is roughly equivalent to the existence 
of a strictly positive SDF. For no arbitrage opportunities to exist, the product of M 
with the value process of any investment strategy must be a martingale []. Under an 
adequate definition of the space of admissible trading strategies, the product MV is 



In a nutshell, an option is an insurance for buying or selling. 

technically, recall that a martingale is a family of random variables £(t) such that the mathematical 
expectation of the increment £(£2) — (for arbitrary t\ < fa), conditioned on the past values £(s) (s < ti), 
is zero. The drift of a martingale is thus zero. This is different from the Markov process, which is better 
known in the physical community, defined by the independence of the next increment on past values. 
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a martingale, where V is the value process of any admissible self-financing trading 
strategy implemented by trading on financial securities. Then, 



V(t) = E, 



M{t) 



(1) 



where s is a future date and Ef [x] denotes the mathematical expectation of x taken at 
time t. In particular, we require that a bank account and zero-coupon discount bonds 
of all maturities satisfy this condition. 

A security is referred to as a (floating-rate) bank account, if it is "locally riskless" 
[]. Thus, the value at time t, of an initial investment of B(0) units in the bank account 
that is continuously reinvested, is given by the following process 

B{t) = £(0)expj^ r{s)ds\ , (2) 

where r(t) is the instantaneous nominal interest rate. 

We further assume that, at any time t, riskless discount bonds of all maturity dates 
s trade in this economy and let P(t, s) denote the time t price of the s maturity bond. 
We require that P(s, s) = 1, that P(t, s) > and that dP(t, s)/ds exists. 

Instantaneous forward rates at time t for all times-to-maturity x > 0, f(t,x), are 
defined by 

/(M) = - a '° g ^ t + l) , (3) 

which is the rate that can be contracted at time t for instantaneously riskless borrowing 
or lending at time t + x. We require that the initial forward curve /(0, x), for all x, be 
continuous. 

Equivalently, from the knowledge of the instantaneous forward rates for all times- 
to-maturity between and time s — t, the price at time t of a bond with maturity s 
can be obtained by 

P(t,s) = exp{- J* f(t,x)dx\ . (4) 

Forward rates thus fully represent the information in the prices of all zero-coupon 
bonds. 

The spot interest rate at time t, r(t), is the instantaneous forward rate at time t 
with time-to-maturity 0, 

r(t) = f(t,0). (5) 

For convenience, we model the dynamics of forward rates. Clearly, we could as 
well model the dynamics of bond prices directly, or even the dynamics of the yields to 
maturity of the zero-coupon bonds. We use forward rates with fixed time-to-maturity 
rather than fixed maturity date. The model of HJM starts from processes for forward 
rates with a fixed maturity date. This is different from what we do. If we use a "hat" 
to denote the forward rates modeled by HJM, 

f(t,s) = f(t,s-t) (6) 



7 A security is "locally riskless" if, over an instantaneous time interval, its value varies deterministically. 
It may still be random, but there is no Brownian term in its dynamics. 
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or, equivalently, 

f(t,x) = f(t,t + x) (7) 
for fixed s. Musiela (1993) and Brace and Musiela fl(| define forward rates in the same 



fashion. Miltersen, Sandmann and Sondermann [171, and Brace, Gatarek and Musiela 



1 18] use definitions of forward rates similar to ours, albeit for non-instantaneous forward 
rates. Modelling forward rates with fixed time-to-maturity is more natural for thinking 
of the dynamics of the entire forward curve as the shape of a string evolving in time. 
In contrast, in HJM, forward rate processes disappear as time reaches their maturities. 
Note, however, that we still impose the martingale condition on bonds with fixed 
maturity date, since these are the financial instruments that are actually traded. 



3 Stochastic strings as solutions of SPDE's 

In a nutshell, the contribution of this paper consists in modelling the dynamical evo- 
lution of the forward rate curve by stochastic partial differential equations (SPDE's) 
||. In the context of continuous-time finance, this is the most natural and general 
extension that can be performed ^|. 

Financial and economic time series are often described to a first degree of approxi- 
mation as random walks, following the precursory work of Bachelier Q and Samuelson 
iflOfl . A random walk is the mathematical translation of the trajectory followed by 
a particle subjected to random velocity variations. The analogous physical system 
described by SPDE's is a stochastic string. The length along the string is the time- 
to-maturity and the string configuration (its transverse deformation) gives the value 
of the forward rate f(t,x) at a given time for each time-to-maturity x. The set of 
admissible dynamics of the configuration of the string as a function of time depends on 
the structure of the SPDE. Let us for the time being restrict our attention to SPDE's 
in which the highest derivative is second order. This second order derivative has a 
simple physical interpretation : the string is subjected to a tension, like a piano chord, 
that tends to bring it back to zero transverse deformation. This tension forces the 
"coupling" among different times-to-maturity so that the forward rate curve is at least 
continuous. In principle, the most general formulation would consider SPDE's with 
terms of arbitrary derivative orders.^ However, it is easy to show that the tension term 
is the dominating restoring force, when present, for deformations of the string (forward 
rate curve) at long "wavelengths", i.e. for slow variations along the time-to-maturity 
axis. Second order SPDE's are thus generic in the sense of a systematic expansion 

In the framework of second order SPDE's, we consider hyperbolic, parabolic and 
elliptic SPDE's, to characterize the dynamics of the string along two directions : inertia 
or mass, and viscosity or subjection to drag forces. A string that has "inertia" or, 
equivalently, "mass" per unit length, along with the tension that keeps it continuous, 

8 Furthcr extensions will include fractional differential equations and integro-differential equations, includ- 
ing jump processes. 

9 Higher order derivatives also have an intuitive physical interpretation. For instance, going up to fourth 
order derivatives in the SPDE correspond to the dynamics of a beam, which has bending elastic modulus 
tending to restore the beam back to zero deformation, even in absence of tension. 

10 There are situations where the tension can be made to vanish (for instance in the presence of a rotational 
symmetry) and then the leading term in the SPDE becomes the fourth order "beam" term. 
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is characterized by the class of hyperbolic SPDE's. For these SPDE's, the highest order 
derivative in time has the same order as the highest order derivative in distance along 
the string (time-to-maturity). As a consequence, hyperbolic SPDE's present wave-like 
solutions, that can propagate as pulses with a "velocity". In this class, we find the 
so-called "Brownian sheet" which is the direct generalization of Brownian motion to 
higher dimensions, that preserves continuity in time-to-maturity. The Brownian sheet 
is the surface spanned by the string configurations as time goes on. The Brownian 
sheet is however non-homogeneous in time-to-maturity, which led us to examine other 
processes. 

If the string has no inertia]^] its dynamics are characterized by parabolic SPDE's. 
These stochastic processes lead to smoother diffusion of shocks through time, along 
time-to-maturity. 

Finally, we mention the third class of SPDE's of second-order, namely elliptic partial 
differential equations. Elliptic SPDE's give processes that are differentiable both in x 
and t. Therefore, in the strict limit of continuous trading, these stochastic processes 
correspond to locally riskless interest rates. 

For the sake of completeness and clarity, we briefly summarize useful facts about 
PDE's (See for instance Ref.[^(| and, in particular, their classification and the intuitive 
meaning behind it. We restrict our discussion to two-dimensional examples. Their 
general form reads 

, (M) ^ +2B(t ,^ +c(t ,^. F(M , /(t , x) ,^,^, s) , 

(8) 

where f(t,x) is the forward rate curve. S(t,x) is the "source" term that will be 
generally taken to be Gaussian white noise r](t, x) characterized by the covariance 

Cov [rf(t, x), 7](t', x')] = 5{t - t') 5(x - x') , (9) 

where 5 denotes the Dirac distribution. Expression (|8[) is the most general second-order 
SPDE in two variables. For arbitrary non-linear terms in F, the existence of solutions 
is not warranted and a case by case study must be performed. For the cases where F 
is linear, the solution f(t,x) exists and its uniqueness is warranted once "boundary" 
conditions are given, such as, for instance, the initial value of the function /(0, x) as 
well as any constraints on the particular form of equation (|8|). 

Equation (||) is defined by its characteristics, which are curves in the (t, x) plane 
that come in two families of equation : 



Adt ={B + VB 2 - AC)alx , (10) 



Adt = (B- VB 2 - AC)dx . (11) 

These characteristics are the geometrical loci of the propagation of the boundary con- 
ditions. 

Three cases must be considered. 
• When B 2 > AC, the characteristics are real curves and the corresponding SPDE's 
are called "hyperbolic". For such hyperbolic SPDE's, the natural coordinate system is 



11 Or if the inertia term is negligible compared to the drag term proportional to the first time derivative 
(so-called overdamped dynamics). 
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formed from the two families of characteristics. Expressing (||) in terms of these two 
natural coordinates A and fx, we get the "normal form" of hyperbolic SPDE's : 

-^L = P(X, (j.) ^ + Q(A, m) ^ + R(\ V)f + S(X, u) . (12) 

The special case P = Q = R = with S(X, fi) = r/(A, //) corresponds to the so-called 
Brownian sheet, well studied in the mathematical literature as the 2D continuous 
generalization of the Brownian motion. 

• When B 2 = AC, there is only one family of characteristics, of equation 

Adt = Bdx . (13) 

Expressing (^|) in terms of the natural characteristic coordinate A and keeping x, we 
get the "normal form" of parabolic SPDE's : 

= K(X, fj) |£ + L(A, ^ + M(X, fj)f + S(X, fj) . (14) 

The diffusion equation, well-known to be associated to the Black-Scholes option pricing 
model, is of this type. The main difference with the hyperbolic equations is that it is no 
more invariant with respect to time-reversal t — * —t. Intuitively, this is due to the fact 
that the diffusion equation is not conservative, the information content (negentropy) 
continually decreases as time goes on. 

• When B 2 < AC, the characteristics are not real curves and the corresponding SPDE's 
are called "elliptic". The equations for the characteristics are complex conjugates of 
each other and we can get the "normal form" of elliptic SPDE's by using the real and 
imaginary parts of these complex coordinates z = u±iv : 

d 2 f d 2 f df df 

7T2+7A = Tjr + U^f + Vf + S. (15) 
ou z ov z ou ov 

There is a deep connection between the solution of elliptic SPDE's and analytic func- 
tions of complex variables. 

We have shown [fl2} | that hyperbolic and parabolic SPDE's provide processes re- 
ducing locally to standard Brownian motion at fixed time-to-maturity, while elliptic 
SPDE's give locally riskless time evolutions. Basically, this stems from the fact that the 
"normal forms" of second-order hyperbolic and parabolic SPDE's involve a first-order 
derivative in time, thus ensuring that the stochastic processes are locally Brownian in 
time. In contrast, the "normal form" of second-order elliptic SPDE's involve a second- 
order derivative with respect to time, which is the cause for the differentiability of the 
process with respect to time. Any higher order SPDE will be Brownian-like in time if 
it remains of order one in its time derivatives (and higher-order in the derivatives with 
respect to x). 

4 No-arbitrage condition : derivation of the gen- 
eral condition 

We now proceed to derive the general condition that the forward rate equation must 
obey to be compatible with the no-arbitrage constraint. 
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From (Q), we get 

d t log P(t, a) = f(t, x) dt - f dy dtf(t, y) , (16) 

Jo 

where x = s — t. We need the expression of -p7^y which is obtained from (^) using 
Ito's calculus In order to get Ito's term in the drift, recall that it results from the 
fact that, if / is stochastic, then 

■w = T** + 2J dx f dx ' mt, W) Cov ld ' m xiMlt - x ' )] • (17) 

where Cov [dtf(t, x), dtf(t, x')] is the covariance of the time increments of f(t,x). 
Using this Ito's calculus, we obtain 



dPjt, s) 
P(t,s) 



dt f(t,x) 



x ~y rx px 

dy E t dt f(t,y) + - / d V dy' Cov [d t f(t,y)d t f(t,y') 
* Jo Jo 



dy[d t f(t,y)-E tAf (t,y)}. (18) 



We have explicitely taken into account the fact that dtf(t, x) may have in general a 
non-zero drift, i.e. its expectation 

E t>dtf (t,x) = E t [dtf(t,x)\f(t,x)} (19) 

conditioned on f(t,x) is non-zero. 

The no-arbitrage condition for buying and holding bonds implies that PM is a 
martingale in time, for any bond price P. Technically this amounts to imposing that 
the drift of PM be zero : 

f(t, x) = f(t, 0) + £ dy EtAf ^ V) ~ \ £ dy £ dy' c(t, y, y') + o(l) , (20) 

assuming that dtf(t, x) is not correlated with the stochastic process driving the pricing 
kernel and using the definitions 

c(t, y, y')dt = Cov[d t f(t, y)d t f(t, y')\ , (21) 

and r(t) = f(t,0) as given by (g). In (pC|), the notation o(l) designs terms of order dt 
taken to a positive power. 

Expression ( p0| ) is the fundamental constraint that a SPDE for f(t,x) must satisfy 
in order to obey the no-arbitrage requirement. As in other formulations, this condition 
relates the drift to the volatility. 

It is useful to parametrize, without loss of generality, 

JWM) = g/M + fc(t)30i (22) 

where hit, x) is a priori arbitrary. The usefulness of this parametrizatio n (p2^ ) stems 
from the fact that it allows us to get rid of the terms f(t, x) and f(t, 0) in fl20|). Indeed, 
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they cancel out with the integral over y of ~^77r~^ • Taking the derivative with respect 
to x of the no-arbitrage condition (p0|), we obtain 

h(t, x) = -- f [c(t, y, x) + c(t, x, y)\ . (23) 



2 

In sum, we have the following constraint that the SPDE for f(t, x) must satisfy 

Et [d t f(t,x)\f(t,x)] = dt 21^.-1 £ dy (cov[d t f(t,y)d t f(t,x)]+Cov[d t f(t,x)d t f(t,y)] 

(24) 

which, for symmetric covariance, lead to 

E t [d t f(t, x)\f(t, x)] = dt - J" dy Cor[dtf(t, y)d t f(t, x)\ . (25) 



This expression (25) is reminiscent of the fluctuation-dissipation theorem in the 



Langevin formulation of the Brownian motion [21|, linking the drag coefficient (analog 
to the l.h.s.) to the integral over time of the correlation function of the fluctuation forces 
(analog of the second term of the r.h.s.). In the usual fluctuation-dissipation theorem, 
the drag coefficient is determined self-consistently as a function of the amplitude and 
correlation of the fluctuating force in order to be compatible with the equilibrium 
distribution. Similarly, here the no-arbitrage condition determines self-consistently the 
conditional drift from the covariance of the fluctuations. In contrast, notice however 
that the same quantity / enters in both sides of (|2~4| ) and (|25|). Thus, the no-arbitrage 
condition imposes a condition on the structure of the partial differential equations that 
govern the dynamical evolution of the forward rates f(t,x). 



5 General structure of the SPDE's compatible 
with the no-arbitrage condition 

To derive the structure of the SPDE's compatible with the no-arbitrage condition 



(|24|), we come back to the formulation of [12] and parametrize the time increment of 
the forward rate as 

d t f(t,x) = a(t,x)dt + a(t,x)d t Z(t,x) , (26) 

where Z(t, x) is a infinite dimensional stochastic process which is continuous in x and 
t and a(t, x) and a(t,x) are a priori arbitrary functions of f(t,x). 

Using this formulation (|26|), we have shown previously ]l2| that the condition of no 
arbitrage leads to the following dynamics for the forward rates 

dtf(t, x) = dt + A(t, x)) +a(t, x) d t Z(t, x) , (27) 



where 



A(t,x) = a(t,x)(J dy a{t,y)c z (t,y,xy\ , (28) 
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The first term of the r.h.s. is a drift term that disappear by a galilean transformation of frame. 
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and 

c z (t,y,y') dt = Cov[d t Z(t,y),d t Z(t,y')] . (29) 

We stress that, by construction, the form ( p7[) with (pq) and (|29| ) automatically satisfies 
the condition (|24|). 

In order to get the form of the allowed SPDE's for f(t, x), we recall the requirements 
that are convenient to impose Z, without loss of generality. 

1. Z(t, x) is continuous in x at all times t; 

2. Z(t, x) is continuous in t for all x; 

3. Z(t, x) is a martingale in time t, E [dfZ(t, x)] = 0, for all x; 

4. The variance of the increments, Var [dfZ(t, x)], does not depend on t or x; 

5. The correlation of the increments, Corr [dtZ(t, x), dtZ(t, x')], does not depend on 
t. 

Then, a fairly general stochastic function Z(t, x) obeying these requirements is such 
that (Santa-Clara and Sornette, 1997) 

d t Z{t, x) = dt [ JX dy rj{t, y) , (30) 



where rj(t,x) is a Gaussian white noise characterized by the covariance (||). This leads 
to the following covariance function for Z(t, x) : 



cz(t,y,y') = Jj^ for j(x)<j(x') , (31) 



and the roles of x and x' in (^Tj) are reversed if j(x) > j(x'). 
Inserting (|30|) in (|27|) , we get 



df(t,x) 8f(t,x) f [* j(y)\ a(t,x) fiW 

x /j(x) 

Multiplying both sides of this equation by ^ ~( and taking the partial derivative with 
respect to x finally yields 



d_ 

dx 



VJjx)f df(t,x) df(t,x) \] d r r— [* j(y) 
^x)[—dt dx-)\ = dx-[v j{x) J dya{t ' y) VJ(x) 



3{x) 



(33) 

where we have used that r](t,j(x)) = rj(t, x) /\/\dj(x) jdx\ . This provides a first class 
of SPDE, which is in general non-linear since the volatility a(t,x) can be an arbitrary 
function of f(t, x). 

An even more general class of SPDE's for f(t,x) is obtained by using the most 
general stochastic function Z(t, x) obeying the requirements 1 — 5 (Santa-Clara and 
Sornette, 1997) : 



d t z(t,x) = dt -±=J^ x) dy vfov) ■ ( 34 ) 
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The correlation of the increments is 



cz(t,x,x') 




if j(x) < j(x') . 



(35) 



The role of x and x' in (35) are inverted if j(x) > j(x'). This provides a generalization to 
( |30| ) since a different function appears in the correlation function and in the inequality 
condition on x and x' . 

Inserting (|4|) in (|27|), we get 



df(t,x) df(t,x) 



dt 



Ox 



Kv)\ cr(t,x) 



(36) 



Multiplying both sides of this equation by ^-jy and taking the partial derivative with 
respect to x finally yields 



d_ 

dx 



a(t, x) V dt 



df(t,x) 
dx 



d_ 

dx 



\Jj(x) / dy<r(t,y)\ 



l(x) 



! dl(x) idj{x) 



dx 



dx 
(37) 

This provides the most general class of SPDE describing the dynamical evolution of the 
forward rates f(t,x) due to the interplay of a stochastic forcing r) and non-linearities 
embodied in a(t, x). 

According to the classification briefly summarized above, these equations (|33| ) and 
(|3?]) are both of the hyperbolic class. It is noteworthy that the no-arbitrage condition 
excludes the parabolic class. Physical strings obey hyperbolic equations for zero or 
small dissipation and parabolic equations in the over-damped limit. A posteriori, it is 
not surprising that we find that the no-arbitrage condition, which implies the absence 
of market friction, corresponds to the first class. Notice that the term d/foaO 

can be replaced by 9 ^g^ 



rj(t,x) 



dt 



9f(t,x) 
dx 



l.h.s. of (|3 



become of the form 



when changing the variables (i, x) to (t,s = t + x) and the 

VW) df(t, s ) 

a(t,s) dt 

These equations (^) and ( |3~7|) are characterized by partial derivatives and two 
source terms in the r.h.s., the first one being locally adapted and the second one 
corresponding to the influence of external noise. These equations can be explicitely 
solved if a is specified and is independent of /(t, x), thus keeping the equations linear, 
as done in (Santa-Clara and Sornette, 1997). The problem of finding a solution of 
these equations when a(t,x) depends on f(t,x) is much more complex and belongs 
to the vast class of generally unsolved non-linear partial differential equations. Such 
equations have been encountered in many different areas, in particular some apparently 
particularly simple nonlinear PDE's have been found to exhibit the most complex 
phenomenology one can imagine. A vivid example is the Navier-Stokes equation of 
fluid motion leading to fluid turbulence when its single parameter (viscosity) goes to 
zero |2^| . Extrapolating on these superficial analogies, we conjecture that it might be 
possible to postulate a simple form for the nonlinear dependence of a as a function of 
f(t, x), such that one or two real parameters might embody the full phenomenology of 
observed forward rate statistics. 



11 



6 A parametric example 



Empirical observation shows that correlation between two forward rates, with maturi- 
ties separated by a given interval, increases with maturity. We have shown previously 
fl~2|| that the parametrization j(x) = e^ x ' leads to the simple and general condition 
that i{x) must be a function which increases more slowly than x, i.e. i(x) must be 
concave (downward). In order words, j{x) must grow slower than an exponential. Any 
expression like j(x) = exp(kappax a lpha) with an exponent < a < 1 qualifies. This 
corresponds to 

c z (t,x,y) = e-^ a -y^ . (38) 
Reporting this choice in (|33|) yields the following SPDE : 

e {K/2)xa f df{t,x) _ df(t,x) Y = 
a(t,x) \ dt dx J _ 

+V^x^e {K/2)xa rj(t,x) . (39) 

7 Reduction to TV-factors models 

Let us now show how our model naturally encompasses the usual HJM formulation of 
forward interest rates in terms of N factors. The idea is that N factors models can be 
obtained as truncations at the order N in a way similar to a galerkin approximation, 
corresponding to a finite "resolution" , of an infinite expansion over the eigenfunctions 
of the operator defining the partial differential equation. This point of view allows one 
to clearly see both the relationship with previous formulations and their limitations. 
It also gives a justification to N factor models in the following sense. As the stochastic 
string description constitutes the most general description of the forward rate curve, 
the fact that the N factors models emerge naturally by a truncation of this general 
formulation, amounting to limit the resolution of the fluctuations in the time of matu- 
rity axis, justifies their mathematical status as simply a degree of approximation of an 
ideal general description. 

To keep the exposition as general as possible, let us call C the linear operator defined 

by 

Cf(t,x) = V (t,x) . (40) 

We simplify the problem by restricting our analysis to linear SPDE's. This allows us to 
use the general property of linear operators that the Green function G(t, x\v, y) defined 
by 

CG{t, x\v, y) = 5(t - v) 5(x - y) (41) 

can be expressed as |20|l 

oo 

G(t, x\v, y) = J2 *) Mv) ■ ( 42 ) 

n=l 

The n(x) are the eigenfunctions of C The <j) n {x) are the eigenfunctions of the adjoint 
operator £*, which is the same as C for the hyperbolic string wave equations considered 
above (self-adjoint case) (but would not be the same for a parabolic string equation). 



d_ 

dx 



d 

dx 



dy a(t,y)e iK/2)y ° 
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The l n (v,t) are a set of functions which depend on C Now, the eigenfunctions of C 
form an orthogonal basis on which one can expand the source term r](t, x) : 

oo 

v (t,x) = 1>j(x) ■ (43) 

3=1 

The delta-covariance property of rj(t, x) allows us to choose that of the T)j(t) as follows : 

Cov[r lj (t),r lk (t')]=g jk (t)5(t-t') , (44) 

where gjk(t) depends on the form of the operator. Using the general solution of (|4C|) 
which reads 

rt rco 

f(t,x) = f(Q,x) + / dv dy G(t,x\v,y) rj(v,y) , (45) 



o 



we obtain 

rt roo 00 



f(x,t) = dv dy^l n {v,t) if) n {x) <pn{y)^Vj( v ) <t>j{y) ■ ( 46 ) 

JO J-oo n=1 j=1 

By the orthogonality condition J^° oc dy(j) n (y) ipj(y) = S n j, we obtain 

f(x,t) = / dvY^l n (v,t) Mx)r]n(v) . (47) 

We recognize a iV-factor model by truncating this sum at n = N and with the identi- 
fication 

\Jgnn{t) l n (v,t) lp n (x) = <J n (t,x) , (48) 

and 

^ dt ~ d t W n (t) . (49) 



V gnn{t) 

Since the functions i^ n (x) are more and more complicated or convoluted as the order 
n increases, the truncation at a finite order N indeed corresponds to a finite resolution 
in the driving of the forward rate curve. For instance, ipi(x) can be a constant, ip2(x) 
can be a parabola with a single maximum, ips(x) can be a quartic (two up maxima 
and one down maximum), etc. 

8 Option Pricing and Replication 

In this section, we derive the general equation for the pricing and hedging of interest 
rate derivatives. Our derivation is restricted to linear SPDE's, i.e. to the cases where 
a(t,x) is not an explicit function of the forward rates /. 

8.1 The general case 

In general, European term structure derivatives in our model have a payoff function of 
the form 

C(t,{P(s)},r(t)) , (50) 
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where the argument {P(s)} indicates that the payoff depends on the price at time s 
of all bonds of different time of maturities, i.e. on the full forward rate curve -P| The 
problem we address is that of hedging this claim by trading bonds. Hedging these 
claims in general implies trading in an infinity of bonds, so that the claim's price at 
time t will be a function of the entire forward rate curve. 

There is thus in general no hope of being able to perfectly hedge interest rate 
contingent claims with a finite number of bonds. We therefore extend the space of 
admissible trading strategies to include density valued portofolios. 

We denote by h(t, s) the density (number) of bonds held in the portfolio at time t, 
with maturity s, and let g(t) be the amount invested in the bank account. Then, the 
value of an investment strategy is 

/•oo 

V(t) = g(t)B(t) + J h(t,u)P(t,u)du 

Let us call V\ the value of the portfolio made of the contingent claim, with value 
process C, financed at the risk-free rate r(t), and V2 that of the replicating portfolio 
consisting of bonds and of money invested in the bank account at the risk-free interest 
rate r(t) (equal to the the spot rate f(t, 0) which may fluctuate but is risk-less in the 
sense that it gives the instantaneous payoff of cash invested in the bank account). At 
time < t < T, the variation of V\ with time, discounted by the risk-free interest rate, 
is 

d t V 1 (t) = d t C(t) - r{t)C{t)dt . (51) 
The variation of V2 discounted by the risk-free interest rate is 

/■oo 

d t V 2 (t)= du h{t,u)[d t P(t,u) -r(t)P(t,u)dt] . (52) 



There are no other terms in ( p2| ) as a sale or purchase of bonds and deposit or with- 
drawal on the bank account correspond only to a change of the nature of the investment 
but not to a change of wealth. The maturity of a bond is not a change of wealth either. 
Only the variations of the bond prices have to be taken into account. Note that the 
term dfB(t) — r(t)B(t)dt vanishes identically from the definition of the spot rate and 
thus the amount invested in the bank account does not contribute to the discounted 



variation of V2. The standard replication argument [23| for option pricing and hedging 
amounts to equate the variations in values of the two portfolios. We need some more 
ingredients before carrying out this program (see the appendix for a pedagogical expo- 
sition of this formulation for the standard European call option problem) . See Ref . p^ j 
for a different approach in the non-Gaussian case. 

We assume that bond prices are driven by one of the stochastic string processes 



Z(t,x) introduced in [12] and used above. C(t) is a priori a function of a continuous 
infinity of bond prices at time t. The relevant mathematical tool to calculate dtC is 
that of functional derivation. We have, up to order dt, 



, _ dC , f°° , dC , „, , 
d t C = —dt + / du —-—d t P(t, u) 
dt Jt oP(t, u) 



dP(t, u) 

13 Expressing the payoff as a function of the full forward curve is the same as expressing it as a function of 
the corresponding continuous set of bonds. 
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I roo roo d^C 

+ U dU l dv d P(t,u)dP(t,v) C ^ PM ^ P ^ 
dC 1 d 2 C 1 f°° d 2 C 

+ w dtr + 2^ YaT[dtr] + 2j t du iMm^f m[dtP{t ' u) ' dtr{t)] ■ (53) 

In we have taken into account that the option price C is also a function of the 
stochastic spot rate r(t). 
d t P(t, s) is given by @ 

r rs-t i 

(54) 



d t P(t,s) = dt P(t,s) 



f(t,0)-P(t,s) f dya(t,y) d t Z(t,y) 
Jo 

This expression can be derived directly from (|lG|), (p7|) and Ito's calculus. Thus 

Cov[d t P(t,n),d 4 P(t,w)] = 

P(t,u) P{t,v) f ' dya(t,y) f * dy' ff (t, y')Cov [d t Z(i, y), d t Z(i,y')] • (55) 
jo Jo 

We also have 

Var[d t r] = [a(t, 0)] 2 Var[d 4 Z(t, 0)] , (56) 

and 

Cov[tkP(t,*),dtr(t)] = -(r(t,0) P(t,s) T *dy a(t,y) Cov [d t Z(t, y), d t Z(t,0)} . 

Jo 

(57) 

The portfolio of bonds replicates the option if dVi(t) = dV^it). We already see that 
a necessary condition for the replication to be perfect (the market to be complete) is 
that the replicating portfolio is made of a continuous infinity of bonds. Technically, 
this comes from the fact that the time differential of the option price C leads, by the 
functional differential, to a continuous integral over all bonds with time-to-maturity 
larger than t when the payoff is indeed dependent on the continuous infinity of bond 
maturities. 

The stochastic part proportional to dtP(t,u) in the replicating equation dV\(t) = 
dViif) cancels out if we choose 

h ^" ) = wBr)- (58) 

which is the usual delta hedging. But this is not enough as the stochastic term pro- 
portional to dtr still remains. Prom the fact that r(t) = f(t, 0) and by definition 

P(M) = exp{- T 1 dy f(t,y)} , (59) 
Jo 

we see that bonds close to maturity are driven by the same stochastic innovations as 
the spot rate. To make sense of this statement, one has to be careful on how the limit 
s — > t is taken. The standard way to tackle this problem is to remember that the 
continuous formulation is nothing but a limit St — > of discrete time increments. We 
thus have, instead of (|59"|), 

P(t, s) = exp{-5t [f(t, 0) + f(t, St) + f(t, 2St) + ... + f(t, (n - l)St)]} , (60) 
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where n = (s — t)/5t. In particular, 



p(t, t + 5t) = exp{-«5t f(t, 0)} = 1 - St f(t, 0) , (61) 

where the second equality becomes asymptotically exact for very small St. This shows 
that it is possible in principle to hedge the spot rate by bonds that are very close to 
maturing. From (|6l|), we obtain 

d t P(t, t + St) = -St d t f(t, 0) = -St d t r(t) . (62) 

From this, we see that the stochastic part proportional to dtr(t) in the replicating 
equation dV\{t) = dV^it) cancels out if we add the quantity of bonds 

Shit - t+dt '> = -iW' <63) 

to the previous quantity h(t,t + St) = gp^j^g^ obtained from (|58|) for the bonds going 
to maturity. Notice that the two quantities are equal since one can replace —5tdr(t) 
by dP(t, t + St) as seen from (|62[). 

To summarize, the hedging strategy is given by 

h(t,u) = \2-Y(u-St)]^L, (64) 

where Y(x) is the Heaviside function equal to 1 for x > 1 and zero otherwise. Thus, we 
recover the usual delta hedging, except for the factor 2 for bonds close to maturation 
which results from the existence of a stochastic spot interest rate. 
The deterministic equation of the option price is then 



poo roo d C(t) 1 f°° a C 

I du l dvA ^ v) ^,^P(M) P( ^ )P( ^ )+ 2l ^^PM^^'^'^ 
Yai[d t r] d 2 C dC(t) . , , [°° ro dC(t) 



+ —^-r(t)C(t)+r(t) / du [2-Y(u-St)] / \ P(t,u) = (65) 



2 dr 2 dt w w w J t 1 v n 8P(t, u) 
where 

ru—t rv—t 

A(t,u,v)= dy a(t,y) dy a(t,y')c z (t,y,y') . (66) 
J o J o 

This equation is correct only when <r(t, is not an explicit function of the forward 
rates /, a case corresponding to linear SPDE's (|37|). When this is not the case, i.e. 
when a(t, x) is a function of /, we see from (|65| ) with (|66|) that C should also be a 
function of /, in addition to be dependent on t,{P(s)} and r(t). As a consequence, 
the total time derivative dtC must contains terms involving partial derivatives with 
respect to / that must be treated self-consistently with the other terms in (|53|). 



8.2 Bond derivatives 

The general case simplifies greatly when the contingent claim has a payoff that can be 
written as a function of a single bond price. We see that pricing and hedging these 
derivatives is much easier. This is very interesting since most derivatives that are 
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actually traded have payoffs that can be written as a function of a finite number of 
bond prices. 

The simplest case of a bond option is a European call or put on a zero-coupon 
bond. Consider a call with maturity s on a bond with maturity s + r, with strike price 
K. The payoff at maturity is the amount 

C(s, P(s, s + t))= max(P(s, s + r) - K, 0) , (67) 

The price of this claim at time t < s will only be a function of P(t,s + r), and the 
derivative can be hedged with a position in this bond alone. 

In the same manner, caps [^], floors 0, collars ™ and swaptions q] have payoffs 
that, in general, can be written as a function of the prices of a finite set of bonds. The 
simplest representative example of this type of options is a caplet. A caplet (settled 
in arrears) pays at maturity, the maximum of the LIBOR (London Interbank Offered 
Rate) rate for the caplet period (set at the beginning of the period) minus the cap rate 
and zero, multiplied by the principal amount. Let the current time be t, the beginning 
of the caplet period be s and the length of the caplet period be r. Denote the principal 
amount by V, let the LIBOR rate be L(s,s + t), and the cap rate be K. Then, the 
payoff at date s + r will be 

C(s + t)=Vt max(L - K, 0) 
that we can express in our notation as 

C(s + t) = V max Uo - 1 - Kr, o) 



= V max — r -1-Kt,0 

\P(s,s + t) 

Finally, this payoff is known at time s, so we can write it as 

C{s) = V max (1 - (1 + Kt)P(s, s + r), 0) 

Again, the payoff only depends on the price of a single bond, so that the claim can be 
priced and hedged with that bond. 

If the time t value of the contingent claim depends only on t and P(t,s), the 
replicating portfolio can have only the bond P(t,s). Then, the previous calculation 
simplifies as all functional derivations transform into the usual derivation and we get, 
instead of (^) : 

1., , rr ,, , l2 d 2 C 1 d 2 C „ , , . \w[d t r]d 2 C 

r)C 8C 8(1 

+ r(t)p(t ' s) 8pm + r(f)p(t ' 1 + 5t) mtT^W> + * ~ rif)c = (68) 



1 A cap guarantees a maximum interest rate for borrowing over a determined time horizon. 

15 A floor guarantees a minimum interest rate for an investment over a determined time horizon. 

16 A collar is a contract in which the buyer is guaranteed an interval of interest rates, with a maximum 
rate. Its sale is thus the association of the buy of a cap and the sell of a floor. 

17 A swaption is an option on a swap. Simply put, an interest rate swap is a contract between two parties 
that exchange for a determined time two interest rates, for instance a short-term and a long-term. 
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where 

PS — t P 8' — 't 

A(t,s)= dy<r(t,y) dy' a(t,y') c z (t,y,y') . (69) 
Jo Jo 

The equation ( |6^ ) is similar to the usual Black-Scholes equation but has two differences. 
First it has a time-dependent diffusion coefficient A(t, x). But more importantly, it 
shows that the option price is function of two bond prices, the bond underlying the 
writing of the option and the bond just before maturation. We see that the expression 



(65) has a priori no mathematical sense in the continuous limit but takes full sense 
when we replace all derivatives by their discrete differences. This is a novel situation 
brought about by the structure of our model defined in terms of correlated but different 

shocks for each maturities as seen from equation d27j). In fact, the continuous limit can 

dC 
dP(t,t+St) 

dC 



be retrieved by remarking that ap/? 1 ?, ^ should be of order St : 



dP(t,t + 5t) 

8C 
3 (t,t+<5t) 



St w(t,P{t,s)) . (70) 

Then, the term r(t)P(t,t + St) 9P ^ +St } in (||) becomes -r(t) w(t,P(t,s)) P(t,t + 
5t) . This situation is similar to a boundary layer for singular perturbation prob- 
lems and in hydrodynamics where a special treatment has to be developed close to a 
boundary, here the spot maturity. 

This equation contrasts with the pricing PDE that is usually presented in term 
structure models. In models with state variables such as CIR, the objective is to price 
derivatives as a function of those state variables and time, and so the PDE is set 
with respect to them. The solution of the PDE we present is not a function of state 
variables, which do not exist in our model, but rather is a function of the price of the 
bond underlying the derivative. 



It is interesting to compare this pricing equation (68) with (69) to the corresponding 
pricing equation for the standard HJM model with a single Brownian motion driving 
the full forward rate curve. In this case, we have 



dtP(t,s)=P(t,s) 

Thus 



s-t 

/ a(t,y)dy 
o 



dt + dW(t) 1° dy a(t,y) . (71) 
o 



Var [d t P(t, s)} = dt [P(t, s)} 2 ^ S * a(t, y) dy^j ' . (72) 

Using the same method as above, except for the term r(t)P(t,t + St) gpn^m which 
is absent, we get the same pricing equation (|6§| ) with the instantaneous "diffusion 
coefficient" 2 

A(t,s)=([ S a(t,y)dy) , (73) 



o 



instead of (|69"D. The expression ( |69"D reduces to ( |73[ ) for cz = 1 as it should, i.e. 
for perfect correlations along the time-to-maturity axis, which corresponds to a single 
Brownian process (single factor) driving the whole forward rate curve. 
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9 Conclusion 



All previous models of interest rates envision the fluctuations of the forward rate curve 
as driven by one or several (multi-factors) "zero-dimensional" random walk processes 
acting either on the left-end-point of the string, the short-rate, or on the whole function 
f(t,x) simultaneously. To sum up the situation pictorially, in the first class of models, 
the forward rate curve is like a whip held by a shaking hand (the central banks!?), 
while in the second class, imagine a large hammer of the size of the full curve hit- 
ting it simultaneously through an irregular pillow (the volatility curve) transmitting 
inhomogeneously the impact along the curve. 

The present paper, which extends Ref.[^], explores the more general situation 
where all the points of the forward curve are simultaneously driven by random shocks. 
Pictorially, this is like a string submitted to the incessant impacts of rain drops. We 
have argued that the condition of absence of arbitrage opportunity provides a con- 
structive principle for the establishment of a general theory of interest rate dynamics. 
We have derived the general condition (P5|) that a stochastic partial differential equa- 
tion (SPDE) for the forward rate must obey. Using our previous parametrization in 
terms of string shocks |12| , we have derived the general structure ( |37| ) of the SPDE for 
forward rates under the no-arbitrage condition. It is noteworthy that the no-arbitrage 
condition excludes the parabolic class of partial differential equations, i.e. those that 
usually describe the physical strings submitted to strong fluctuations! 

In a second part, we have derived a general approach to price and hedge options 
defined on bonds and thus on the forward rate curve, in terms of functional partial 
differential equations. We have found that the usual Black-Scholes replication strategy 
can be adapted to this situation, provided that special terms be added to account self- 
consistently for the instantaneous interest rate. Our approach is limited to the case 
where the volatility <j(t, x) does not depend on fit, x) itself, thus excluding the a priori 
most interesting class of stochastic non-linear partial differential equations. 

Let us end with a conjecture. It may be conceivable that a simple form for the 
nonlinear dependence of a as a function of f(t,x), with one or two real parameters, 
might embody the full phenomenology of observed forward rate statistics, hence con- 
stituting a truly fundamental theory of forward rate curves. It would be fundamental 
in the sense that the properties of the shocks would be generated dynamically by the 
nonlinear interactions similarly for instance to the coherent vortices in turbulence for 
instance, and in contrast to the usual external shock assignments in the existing linear 
models. The complex behavior would then result from the interplay between external 
factors represented by the noise source and the non- linear dynamics. Persuing along 
this conjectural tone, this approach might provide an line of attack for understanding 
the recent empirical finding p5| of a causal information cascade across scales in volatil- 
ities, occurring from large time scales to short times scales in a way very similar to the 
celebrated Kolmogorov energy cascade proposed for fluid turbulence. 
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Clara for helpful discussions and to D. Stauffer for a careful reading of the manuscript. 
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APPENDIX 

The simplest option pricing problem (the so called 'European call options') is the 
following : suppose that an operator wants to buy from a bank a given share, a certain 
time t = T from now (t = 0), at a fixed 'striking' price x c - If the share value at 
t = T, x(T), exceeds x c , the operator 'exercises' his option. His gain, when reselling 
immediately at the current price x(T), is thus the difference x{T) — x c . On the contrary, 
if x{T) < x c the operator does not buy the share. What is the price C of this option, 
and what trading strategy <f> should be followed by the bank between now and T, 
depending on what the share value x(t) actually does between t = and t = T? 



In the standard treatment in continuous time [23], one forms the portfolio F = 
—C + x(j) such that dtF remains indentically zero. Here, we propose a slightly different 
derivation of Black- Scholes' result based on the idea that the bank constructs a portfolio 
that replicates the option exactly, thereby eliminating all risks. This is of course only 
valid under the restricted assumptions of continuous trading, Gaussian random walk 
of market prices, and absence of market imperfections and transaction costs. 

The idea is to compare the two points of view of the option buyer and of the option 
seller. During the time increment dt, their respective change of wealth discounted by 
the risk-free interest rate r is 

dW a = dC- rCdt , (74) 

for the buyer who owns only the option and 

dW v = 4>[dx - rxdt] , (75) 

for the seller who possesses a portfolio made of (ft underlying stock shares. dC — rCdt 
is the gain or loss of the buyer above the risk-free return. cf)[dx — rxdt] is the gain or 
loss of the seller for a price variation of the stock above the risk-free return. The fair 
price of the option and the hedging strategy that the seller must follow are those such 
that the return is the same for both traders, namely 

dW a = dW v . (76) 

This equality makes concrete the fundamental idea of Black and Scholes that the fair 
price and the hedging strategy are univoquely determined for all traders, independently 
of their risk aversion, if the seller can replicate the option. 
Using Ito's formula to calculate dC, we have : 

dC(x,x c ,T-t) ^ dC(x,x c ,T-t) j D d 2 C(x,x c ,T-t) ^ 
^ = m dt + Wx dx + — 2 dt. (77) 



Inserting this expression in (76) shows immediately that, if one makes the choice 

, ,* dC(x,x c ,T -t) 

<P = <t> = 5 , 

OX 



(78) 



then the only stochastic term, namely dx, cancels out! The comparison between the 
two returns of the buyer and seller become certain. In this case, the equation (f76|) 
provides the following deterministic partial differential equation for C: 

dC(x,x c ,T-t) dC(x,x c ,T-t) D d 2 C(x,x c ,T - 1) „. m . 
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with bouncary conditions the value of the option at maturity C(x,x c ,0) = max(x — 
x c , 0). The solution of this equation is the celebrated Black and Scholes formula JlC|] . 
The hedging strategy is then the derivative of this solution with respect to x. This 
derivation shows very straightforwardly why the average return of the underlying stock 
does not appear; it has been avoided by the replication condition (|76|) with (|78|), 



21 



References 

[1] Brealey, R.A. and S.C. Myers, 1991, Principles of corporate finance, fourth edition 
(McGraw-Hill, New York). 

Bachelier, L., 1900, Theory of speculation (Paris: Gauthier-Villars) 
Mandelbrot, B.B., 1963, Journal of Business 36, 394. 

B. Dubrulle, F. Graner and D. Sornette, eds., Scale invariance and beyond (EDP 
Sciences and Springer, Berlin, 1997). 

Da Prato, G., and J. Zabczyk, 1992, Stochastic equations in infinite dimensions 
(Cambridge University Press; UK). 

Witten E., 1996, "Reflections on the fate of spacetime", Physics Today 49, 24-30. 

Modigliani, F. and M.H. Miller, 1958, "The cost of capital, corporation finance 
and the theory of investment", American Economic Review 48, 655-669. 

Miller, M.H. and F. Modigliani, 1961, "Divident policy, growth and the valuation 
of shares", Journal of Business 34, 411-433. 

Bouchaud, J. -P., N. Sagna, R. Cont, N. El-Karoui and M. Potters, Phenomenology 
of the interest rate curve, preprint cond-mat/9712164 



Black, F. and M. Scholes, 1973, "The pricing of options and corporate liabilities", 
Journal of Political Economy 81, 637-654. 

Taleb, N., 1997, "On the fallacy of transaction costs in option replication", 
preprint. 

Santa-Clara, P. and D. Sornette, 1997, "The Dynamics of the Forward Interest 
Rate Curve with Stochastic String Shocks", submitted to Review of Financial 
Studies ( |http: / / xxx.lanl.gov/abs/cond-mat /9801321| ) 

N. Ikeda et al. eds., 1996, Ito's stochastic calculus and probability theory (Tokyo ; 
New York : Springer). 

Heath, D., R.A. Jarrow and A. Morton, 1992, "Bond Pricing and the Term Struc- 
ture of Interest Rates: A New Methodology for Contingent Claims Valuation", 
Econometrica 60, 77-105. 

Duffie, D., 1996, "Dynamic Asset Pricing Theory", Second Edition, Princeton 
University Press. 

Brace, A. and M.Musiela, 1994, "A Multifactor Gauss Markov Implementation of 
Heath, Jarrow, Morton" , Mathematical Finance 4, 259-283. 

Miltersen, K., K.Sandmann and D.Sondermann, 1997, "Closed Form Solutions for 
Term Structure Derivatives with Log-Normal Interest Rates" , Journal of Finance, 
52, 409-430. 

Brace, A., D.Gatarek and M.Musiela, 1997, "The Market Model of Interest Rate 
Dynamics", Mathematical Finance, 7, 127-155. 

Samuelson, P.A., 1965, "Proof that properly anticipated prices fluctuate ran- 
domly" , Industrial Management Review 6, 41-50. 



22 



[20] Morse, P.M. and H. Feshbach, 1953, Methods in theoretical physics, McGrawHill 
Publishing Company, New York. 

[21] Stratonovich, R. L., 1992, Nonlinear nonequilibrium thermodynamics I : linear and 
nonlinear fluctuation- dissipation theorems (Berlin ; New York : Springer- Verlag) . 

[22] Frisch, U., 1995, Turbulence, the legacy of Kolmogorov (Cambridge University 
Press; UK). 

[23] Merton, R.C., 1990, Continuous-time finance (Blackwell Publishers, Cambridge, 
Massachussetts) . 

[24] Bouchaud, J. -P. and D. Sornette, 1994, "The Black-Scholes option pricing prob- 
lem in mathematical finance : Generalization and extensions for a large class of 
stochastic processes", J.Phys.I France 4, 863-881. 

[25] Arneodo, A., J.-F. Muzy and D. Sornette, 1998, "Direct" causal cas- 
cade in the stock market, in press in European Physical Journal B 
( http:/ /xxx. laQl.gov/abs/cond-mat/9708012D ■ 



23 



